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We present a theory of the phonon-assisted nonlinear dc transport of 2D electrons in high Landau levels. The 
nonlinear dissipative resistivity displays quantum magneto-oscillations governed by two parameters which are 
proportional to the Hall drift velocity vh of electrons in electric field and the speed of sound s. In the subsonic 
regime, vh < s, the theory quantitatively reproduces the oscillation pattern observed in recent experiments. 
We also find the 7r/2 phase change of oscillations across the sound barrier vh = s. In the supersonic regime, 
Vh > s, the amplitude of oscillations saturates with lowering temperature, while the subsonic region displays 
exponential suppression of the phonon-assisted oscillations with temperature. 
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A high-mobility two-dimensional electron gas (2DEG) in a 
weakly quantizing perpendicular magnetic field B displays a 
rich variety of magneto-oscillation phenomena. These phe- 
nomena are defined by the ratio of energy parameters and the 
Landau level spacing huc, with Wc = eB/mc, or by ratio of 
the spatial length scales and the cyclotron radius R^ — vp/^c, 
where vp (pp = mvp) is the Fermi velocity (momentum). 

Examples of magneto-oscillations include: (i) the 
microwave-induced resistance oscillations (MIRO) yj-lsl] and 
associated zero-resistance states (ZRS) j^]^, controlled by 
the ratio of the the microwave frequency cu and cyclotron fre- 
quency ujc, Cac = i^/wc; ('0 the Hall field-induced resistance 
oscillations (HIRO) |6,7] and associated zero-differential re- 
sistance states (ZdRS) 18[] in the dc electric field E, gov- 
erned by edc = eE{2Rc)/hijJc\ (iii) the phonon-induced re- 
sistance oscillations (PIRO) J 91 41211 . characterized by the ratio 
Eph = w^/wc, where HlUt^ = 2spp is the energy of an acous- 
tic phonon with momentum 2pp and s is the speed of sound. 
The above oscillations stem from commensurability of energy 
change in photon- (MIRO); disorder- (HIRO), and phonon- 
assisted (PIRO) scattering events and the ftwc -periodic oscilla- 
tions of the density of states (DoS). MIRO and HIRO 1 1 , 2,4- 
ISD, as well as their interplay il3l4l5ll . have been studied in 
detail over the past decade. 

In this paper, we present a theory of a phonon-assisted 
transport in strong dc electric fields. In the nonlinear dc- 
response, our results are consistent with the observed evo- 
lution of PIRO in the regime edc < £ph in recent experi- 
ment ilOtl . We predict that the phase change of oscillations 
across the sound barrier edc = Cph, where the Hall drift veloc- 
ity Vh — cE/B equals to the speed of sound. We also analyze 
the temperature dependence of the phonon-assisted transport 
in different regimes. 

Model. For a homogeneous 2DEG in a region of classically 
strong B, when Wc exceeds the transport scattering rate I/tu, 
the dissipative current jx in response to an electric field E 



(pointing along x direction, ip = 0)is 

jxE = 2vq I W^^^'P^^^'{e)de 



dipdip' 
~4^ 



(1) 



Here i^o = m/27r is the DoS at zero magnetic field (hereafter 
we put h — 1) and P^p.^p'{e) denotes the semiclassical scat- 
tering rate of an electron with energy e with the change in its 
momentum direction from (^ to if': 
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In scattering of static disorder, represented by the first line in 
Eq. (|2]l, electron energy e changes by the electrostatic work 
Wip^ip' — ei?i?c(sin</? — sin(^') for displacement Aii of elec- 
tron cyclotron trajectory in electric field, see Fig.[TJ. The sec- 
ond (third) line in Eq. dU describes processes with change 
in electron electrostatic energy Wip,ip' and absorption (emis- 
sion) of a phonon with energy |a;;p_;^/|, where uj^p-^p' = 
2ppssm{ip — (p')/2 for s <^ vp. The Planck's function 
Nui — l/[exp(|aj|/T) — l] accounts forthe thermal occupation 
of phonon modes, and g is the dimensionless electron-phonon 
coupling constant 1161]. 

Factors 7W(e,e') = v{e)9{e')f{e)[l - f{e')] in Eq. © 
contain the Fermi function /(e) = l/{exp[(£ — /i)/r] + 1} 
and the product of the DoS of initial and final states. The DoS 
in the experimentally relevant region of weak magnetic fields. 



< 1, has the form 



i>{e) = v{e)/v() = 1 — 2Acos(27re/ajc), (3) 

where A = cxp(— 7r/wcTq) and Tq is the quantum scattering 
time off disorder 

Using Eqs. ([T]i-([3]i and assuming Trq > 1, we obtain 
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FIG. 2: r^x'^ vs B for several values of current I varying from fiA 
to 50 ^A in 5 fj,A steps. Traces are vertically offset for clarity. Traces 
(A) [10/iA], (B) [15/iA], and (C) [30/iA] correspond approximately 
to the values of I at which 3rd, 2nd, and 1st PIRO peaks disappear 
before evolving into minima at higher currents. Very similar values 
(A) [10/iA], (B) [14^A], and (C) [28/iA] were found in the experi- 
ment (see Fig. 2 of Ref. [10]). 



FIG. 1: Illustration of electron-phonon scattering processes in differ- 
ent regimes: (a) for edc = 0; (b) for edc < tph; (c) for edc = Cph; 
(d) for edc > Eph. Arrows show the change of the total energy of 
electron due to emission or absorption of a 2pir-phonon accompa- 
nied by the shift of the guiding center of the cyclotron orbit by 2Rc 
along or against the dc field. The arrows are directed along the sound 
cone (solid lines), while the strength of dc field determines the tilt 
of Landau levels (maxima of oscillatory DoS Eq. ([3} are marked by 
yellow stripes). Dashed line marks the position of the Fermi level. 
The position of four regions Fig.[T^-d in the eph ~ edc plane is shown 
in Fig.[T^. Fig. [if illustrates the real-space shift AR of the cyclotron 
guiding center and the change in the electron momentum for scatter- 
ing from ipto if'. 



T°^^ in the experimentally relevant limit T ^ Wc and then 
analyze how these rates evolve with lowering T. 
Main results. At high T > W,U!, 



A(w/2T, W/2T) 



we approximate 
1 in Eq. Q and obtain constant F^ ™^ = 
g^T for the smooth part of electron-phonon transport scatter- 
ing rate. The oscillatory part, which represents the effect of 
the Landau quantization, reduces to 
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cos [27r sin Lp^ (eph — edc cos (^+)] 



where po = 7i't'o/(e^^eTtr) and pn = B/{ecNf.) are the 
Drude and Hall resistivities, N^. — mvpiyQ is the electron 
surface density. The effective scattering rate of disorder was 



obtained in |7]: rdis(edc) 



2A'^F"(7redc), with 



^{^) = Yl Jni^)/'''n and Jn{x) Standing for the Bessel func- 



tions; the scattering rate r 



= E^,7^e™('^-'^') is pre- 



sentedby its angular harmonics l/r„ (tq = Tq). 

Here we focus on the phonon-assisted scattering rates 
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where r|,™^(edc, Cph) and rpjf''-'(edc, Cph) denote the smooth 
and oscillatory components of the electron-phonon scattering 
rate, and A{x,y) = S{x)S{x — y)/S{y) is written in terms 
of S{x) = x/sinh x. We first study the behavior of F^ ™' and 



-i(osc 



The behavior of F^ (edc7 Cph) is different in parametric 



regimes, illustrated in Fig. [T^. In the limit of weak electric 
fields, edc ^ 1, the current is linear in the applied dc field and 



rS^^(o,^ph) 



2A^.g^r[Jo(27reph) - J2(27repi0]. (7) 



The oscillations with eph — i^tt/^c in Eq. (|7]i are 
driven by commensurability between the phonon energy 
Wtt = 2pFS and the period of DoS oscillations cUc, see 
Fig. [Di. For eph > 1, Eq. © gives Fj,°if=)(0, eph) = 
{4X^ g^T /TTy'eph) cos(27reph — 7r/4). This regime was studied 
previously by Raichev 1 17]. The difference in the oscillations 
phase of our result and that of Ref. 1,17.] is due to a different 
choice of the photon scattering form-factor in narrow vs. wide 
quantum wells. 

A strong Hall field tilts the Landau levels and changes the 
commensurability conditions for electron-phonon scattering 
as shown in Fig.[T}5-d, resulting in oscillations with combined 
parameters e± = e^c ± Cph- Equivalently, this effect can be 



viewed as a result of the Doppler shift of the phonon modes in 
the frame moving with the Hall velocity vh in y direction. In 
the moving frame, the electric field is absent while the phonon 
speed changes as s — > s—vh cos ip^ in the argument of cosine 
in the second line of Eq. ^. The latter formulation is partic- 
ularly useful in the case of complicated spectrum of phonons. 

can be approximated as 
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This expression can be simplified in three cases when electron 
drifts with subsonic velocity {vh < s or edc < Cph, Fig.[T]3); 
velocity of sound (vh = s or edc = fph, Fig.[T];); and super- 
sonic velocity {vh > s or edc > Cph, Fig.[TJl): 
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Evolution of PIRO in subsonic sector edc < Cph was 
studied experimentally in Ref. [JOJ, where the differential 

resistivity rxx = OE/djx — pjjdjx/dE was measured. 



The PIRO contribution to the differential resistivity 
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can be evaluated using Eqs. (|4|i and (|5]l as 
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In Fig.|2]we plot rxx'^ , calculated according to Eqs. ^ and 



© with 



ps as a function of _B = Bo /eph for several 



values of the total current / = jxW = /oCdc/eph- Here w — 
50 /im is the Hall bar width, Bq = 2ppsmc/e — 7.67 kG 
and Jo — eNeSw — 223 /iA are the reference field and cur- 
rent; values of the parameters s, A^e, and w were taken from 
Ref. llioll . This theoretical result quantitatively reproduces the 
experimentally observed evolution of PIRO with increasing / 
OlOll without fitting parameters. 

The important feature of Eq. (|9]l is the change in the phase 
of oscillations across the sound barrier v^ = s (see Fig. [T];), 
which calls for an experimental verification. In the supersonic 
sector edc > Cph, see Fig. [TJl, the scattering along electric 
field always results in gain in the Landau level index,while 
scattering against electric field reduces the Landau level index. 

In the extreme limit of strong B or soft phonons Cph <C 1 
(see region (e) in Fig.[Tt), Eq. ^ gives 



ri,r^(6dc,o) 

2A2g2T 



Jo{T^£dc) - 'h (Tredc) - '/o(7redc)^2(7redc). 

(10) 



Here, the electron-phonon scattering behaves as a scattering 
off nearly static impurities and the functional form Eq. ( fTOl i 
coincides with that of HIRO, presented above as i^"(7redc), 
for the specific case of isotropic disorder scattering, tq = tq. 

As temperature T decreases, the phonon-assisted magneto- 
oscillations in differential resistivity become suppressed com- 
pletely in the subsonic regime. However, in the supersonic 
regime only e+ oscillations vanish, while amplitude of e_ os- 
cillations saturates, see Figs.|3t-c. 

This behavior can also be inferred from the analytical re- 
sults for Fp^^'^ , Eq. (|5]l, obtained for low temperatures T < 
2pps, but still T > ujc- For edc < Cph, we have 
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^ph'^ is exponentially small for e_ < 0, but as the sound 
barrier is reached, a non-exponential contribution emerges 
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In the supersonic regime, e_ = edc — fph > 0, we obtain 



p(osc), ^ _ 2A g ^c£ph 

-'^ ph ledcEphj — ^3/2 



TT^e 



3/2 
dc 



(13) 



7ersin27re- + ^e^e^^'p"""^'^ sin27re+ 



2£ph(^,/T . 



A similar behavior also occurs with the smooth component 
of the electron-phonon rate T^^ : it vanishes in the subsonic 
regime and saturates in the supersonic regime, see Fig.[3}l. For 
T = 0, the smooth component 
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At finite 



for edc > eph and FJ,™'^ (edc, eph) = otherwise. 

but low temperature, T^^ (edc, eph) remains small in the sub- 
sonic sector and grows continuously across the sound barrier. 
Such contrast in low-temperature behavior of electron- 
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phonon scattering rates T\ ' and T\ ' in the subsonic and 
supersonic regimes can be understood from the structure of 
Eq. dlj. At low temperatures T < 2pFS, the occupation 
number of phonon modes with energy lo ^ 2pFS is expo- 
nentially small. Therefore, the electron scattering off thermal 
phonons becomes negligible. At the same time, the sponta- 
neous emission of phonons depends only on the combination 
/(£m)[l — /(£/)] of electron occupation numbers /(e) at ini- 
tial (Ein) and final (Sf) energies. For emission processes, this 
combination is also exponentially small in equilibrium, but 
may remain finite in a system subject to electric fields. When 






FIG. 3: (a)-(c): Oscillatory part of the phonon-assisted differential 
resistivity r^x'^ in units 2\^g^pD(-^cTa calculated using Eqs. Q, ([5j 
for the case of high [T — 5uuc, panel (a)], intermediate [T = 0.7wc, 
panel (b)], and low [T = 0.25a;c, panel (c)] temperature. Panel 
(d): Smooth part of phonon-assisted scattering rate FpjJ'^ (edc, fiph) 
at zero temperature, measured in units of g^ujc, Eq. J14t . 



the Hall drift exceeds the speed of sound, Sf for electron scat- 
tering with increase in electrostatic energy is larger than Ei 
even for the process with phonon emissions, see Fig. [T}l, and 
the spontaneous phonon emission takes place. 

Conclusions. We developed a theory of phonon-assisted 
magneto-oscillations in strong dc electric fields. In the limit 
of a weak Hall field, edc ^ 1 ^ Cph, or in the limit of quasi- 
elastic phonon scattering, tph <C 1 ^ Cdc, the theory is con- 
sistent with the known results for PIRO (l^ and HIRO |0]. 
In the nonlinear dc-response edc < Cph, the theory repro- 
duces the evolution of PIRO observed in recent experiment, 
cf. Ref. ilOtl and Fig. |2] We find that the phase of oscilla- 
tions changes across the sound barrier at edc = Cph, where the 
Hall drift velocity vh ~ cE/B reaches the speed of sound s, 
which also warrants an experimental verification. In the super- 
sonic regime edc > Cph, the amplitude of oscillations saturates 
at finite value for very low temperatures T <^ uj-^, while in 
the subsonic sector edc < Cph the phonon-assisted cuiTent de- 
cays exponentially with lowering T, as the number of thermal 
phonons decreases. At zero temperature, the smooth part of 
the phonon-assisted nonlinear conductivity also exhibits satu- 
ration to a finite value in the supersonic regime, while vanishes 
in the subsonic regime. 
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